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Abstract. We study soliton solutions to the DKP equation which is defined by the Hirota bilinear 
form, 

r (-4D^Di+i?4+3n2)T-„.T„ = 24r„_iT„+i, 

\ {2Dt + DIt ZD^Dy) T„±i ■ r„ = " 
where tq = 1. The r-functions Tn are given by the pfaffians of certain skew-symmetric matrix. 
We identify one-soliton solution as an element of the Weyl group of D-type, and discuss a general 
structure of the interaction patterns among the solitons. Soliton solutions are characterized by 
4Af X AN skew-symmetric constant matrix which we call the _B-matrices. We then find that one 
can have M-soliton solutions with M being any number from A'^ to 2A'^ — 1 for some of the iN X 4A'^ 
S-matrices having only 2N nonzero entries in the upper triangular part (the number of solitons 
obtained from those _B-matrices was previously expected to be just N). 



1. Introduction 

We study soliton solutions of the DKP equation 00 El , 

r (-4i?,A +Dl + iDl) r„ • T„ = 24r„_iT„+i , 

(1.1) <^ n = l,2..., 
[ {2Dt + DIt 3DxDy) T„±i ■Tn = 

where tq — I, and Dx^Dy and Dt are the usual Hirota derivatives. For each n, the variables 
u = 2(lnT„)a;a: and — Tn±i/Tn define the coupled KP equation (see for example (3), 

(-4Mf + UxxX + GUUx)x + SUyy = 24:{v+V^)xx , 

+ V^xx + SmW^ T K'^xy + r "a '^^) = • 

This set of equations admits a class of particular solutions, called soliton solutions similar to those 
of the KP equation (see for example |^), and in this paper we describe some properties of those 
solutions. 

Remark 1.1. The DKP equation was given as a first member of the DKP hierarchy in In 0, 
the coupled KP equation was introduced as an extension of the KP equation whose solutions are 
given by the pfaffian form. The DKP equation was rediscovered as the Pfaff lattice describing the 
partition function of a skew-symmetric matrix model in ^ [7] , and was also found as an orbit of 
some infinite-dimensional Clifford group action in P). 

The T„ functions are given by the pfaffians of 2n x 2n skew-symmetric matrices Q„ whose entries 
are denoted by Qij for 1 <i < j <2n with Qj i ~ —Qi.j , 

(1.2) T„ = Pf(Q„) = ^ 0-{h,jl,---,in,jn)QtijiQi2,j2---Qir.,Jr. 

1 = ii < . . . < i„ < 2n 
ik < jk, k = 1,.. . ,n 
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Each coefficient cr{ii,ji, ■ ■ ■ ,imjn) gives a sign corresponding to the parity of the permutation, 

1 2 • • • 2n - 1 2n^ 



Sign 



Here the elements Qij satisfy 



— ( 

dtk 



i,j — Q 



i+k,j 



{i,j+k 



fc- 1,2,. 



where ti = x,t2 = y,t3 = t and others are the symmetry parameters. A realization of Qi j is given 

by 0, 

(j-i) 



Q^ 



for i < j 



where 0'-'"'^ — d'^tp/dx'^ (the same for and the functions (p and ^ satisfy the same equations for 

X and t evolutions, 

d(j) d'^(j) dcj) d^4> 
dy dx"^ ' dt dx^ 
For an example of finite dimensional solution, we consider 



(1.3) 



4>{x,y,t) 



M 

E 

m—1 



M 



aniEm{x,y,t), 1p{x,y,t)^ y braEra{x,y,t) 



m—1 



with some constants a^, bm for m = 1, . . . , M. The function E„i{x, y, t) is the exponential function, 

Ejn{x,y,t) := e^"' with 9ra{x,y,t) = p^nX + pl,y + pl,t + 9^^ , 

where Pm and 9^ are arbitrary constants, and throughout this paper we assume the parameters 
p := {pi,p2, ■ ■ ■ ,Pm) to be ordered as, 

(1.4) pi < P2 < ■ ■ ■ < PM ■ 

With those (j) and ip, the elements Qij, I < i < j < 2n become 



Q^ 



(1.5) 



l<fc<i<Af 



/ ^ bkjipkPiY ^{pj 



P'k )Ek, 



= E 

l<fc<i<Af 

where = akk - a;6fc, and i?fc,; := EkEi = exp(6'fc + 9i). 

As a generalization of this form of Qij, we consider an arbitrary M x M skew-symmetric matrix 
B — (bk.i)i<k,i<M- We then note that the 2n x 2n matrix Q„ in the r- function t„ can be expressed 
as 

(1.6) Qn{x, y, t) = £n{x, y, t)B£n{x, y, tf 

where £n is a 2n x M matrix with its transpose f and it is given by 

(El E2 ■ ■ ■ Em 



(2n-l) 
2 



E 



(2n-l) 
M 



In this paper, we discuss a classification problem of several soliton solutions given by 1)1. 6|l in terms 
of the i?- matrix. Note here that M > 2 for n > 1, since the case with M — 2 gives a trivial solution, 
i.e. Ti — Qi,2 = bi 2{p2 ~ pi)EiE2 and u = 2(lnri)2:2; = 0. Also note that one needs M > 2n for 

U = 2{\nTn)xx ^ 0. 

One should also note that if t„+i = (i.e. — 0), the corresponding solution satisfies the 
KP equation. This condition can be of course obtained from the structure of the i?-matrix. For 
example, if we take M = 3, then T2 vanishes identically (the size of the pfaffian is 4 x 4, but the 
independent exponentials are three or less). This implies that t\ with M — 2, gives a solution of the 
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KP equation, and it gives either one KP soliton solution or a resonant Y-shape KP soliton: With a 
3x3 _B-matrix, we have 

n = Ql,2 = fel,2(P2 - Pl)EiE2 + 61,3(^3 - Pl)EiE3 + 62,3(^3 - P2)E2E3 . 

The function u = 2(lnTi)a;a; gives one KP soHton solution if one of bi_j is zero (with others being 
positive), and Y-shape KP soliton if all bij are positive. For example, with 61^2 = 0, we have 
Ti = (&i,3(P3 - Pi)Ei + &2,3(P3 ~ P2)E2)E3 which leads to 

(1.7) w{x,y,t) := ^Inn =^3 + ^{pi + P2) + ^{p2 - pi)tanh^(6'2 - 6*1). 

The asymptotic values of w then take 



pi + p3, for X — > -00 
P2 +P3, for X 00 



Note here that the one soliton exchange the asymptotic values of (1, 3) := P1+P3 and (2, 3) := P2+P3, 
that is, this one soliton permutes the numbers, (1) ^ (2). We then label one soliton of the KP 
equation as an element of the permutation group W, in this case W — S3, the symmetry group of 
order 3. We denote the one soliton solution (|1.7(l by [1 : 2], and call this type of soliton A-soliton 
( "A" stands for type A Lie algebra which is the underlying symmetry algebra for KP equation) . In 
general, we denote one A-soliton by [i : j], if the function w exchanges pi <-> pj with pi < pj. We 
sometime identify [i : j] as an element of the symmetry group. 

A generic solution of the DKP equation is then obtained for M > 4. In particular, we obtain one 
soliton solution in the case M = 4 with the S-matrix having just two nonzero elements in the upper 
triangular part. For example, we consider the _B-matrix having bi^2 = 1, ^3,4 = 1 and all others 
bij — for i < j. Then 

Tl = Qia = {P2 - Pl)EiE2 + (P4 - P3)E3Ei . 

This gives 

9 1^1 1 

(1.8) w{x,y,t) ^Inn = 2 + 2^P3 + Pi - Pi - P2)tanh-(6'34 - 6*12) 

k=l 

with 9ij = 0i + Oj + In \pi — Pj\. For each asymptotic of x — > ±cx), one of the exponential terms in 
the r-function becomes dominant. In the case of (|1.8|l . we have 



w{x,y,t) 



P1+P2, asx— >— 00, 
P3+ Pii as a; ^ 00 . 



Thus, this one D-soliton exchanges (1,2) ^ (3,4) with the values oi w where (i,j) Pi +Pj- This 
indicates a Weyl-action of D-type: For example, tti 2 G may be expressed as 

7i'i^2 • (ii,«2 : H^h) = {i2,H ■ hih)-, H = l,i2 = 3,ii = 2,^2 = 4. 

We denote this D-soliton by [1,2 : 3,4]. Then identifying [1,2 : 3,4] with the permutation of the 
pairs, (1, 2) ^ (3, 4), one has the relation [1, 2 : 3, 4] = [1 : 3] • [2 : 4] = [1 : 4] ■ [2 : 3] where [i : j] is 
the permutation (i) ^ (j) of A-soliton. This relation implies a resonant bifurcation of one D-soliton 
into two A-solitons (see Section 0}. Since one D-soliton has four p,;-parameters, we need to have 
M — AN parameters to describe N D-soliton solution as a solution given by the r-function tn- 
However, we show in this paper that 4A^ parameter solution can contain up to 2A^ — 1 number of 
D-solitons (this is quite different from the solitons of the KP equation). 

In this paper, we study soliton solution of the DKP equation consisting of those A- and D-type 
soliton solutions. In Sectional we give a general structure of the r-functions. Then in Section |21 
we discuss some details of one D-soliton solutions given by the B-matrix of size 4x4. Here we 
also classify the soliton solutions obtained by 4 x 4 i3-matrices. In Sectional we present the soliton 
solutions for the case of 8 x 8 _B-matrices, and classify the soliton solutions consisting of only D-types. 
It turns out that the case with 8x8 i3-matrix can have either two or three D-solitons depending on 
the values of the parameters {pi : i — 1 ...?>} . The generic solution given by the B- matrix having 
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all nonzero entries is then given by four A-solitons. Finally we discuss the general case of AN x AN 
_B-matrices, i.e. M = AN, in Section [S] Then we show that the number of D-solitons for those 
_B-matrices having 2N nonzero entries in the upper triangular part can be any number from N to 
2N — 1, contrary to the previous study (see 0) where the number of D-soliton is expected to be 
just N. 

In jH] , we showed that the A^-soliton solutions of the KP equation can be classified by the Schubert 
decomposition of the Grassmannian Gr{N, 2N). We also expect that iV-soliton solutions of the DKP 
equation can be classified by the similar decomposition of the orthogonal Grassmannain OGr(M, 2M) 
for some M. However, in this paper, we just present elementary feature of the iV-soliton solutions 
based on the pfaffian structure of the r function, which is already complicated but has several 
interesting aspects. We plan to discuss a classification problem based on the geometric structure of 
the orthogonal Grassmannain in a future communication. 

2. Structure of the t-functions 

We first note that the r-function (|1.2|l with the Q„ matrix 1)1. 6|l has the following expansion 
theorem, i.e. the pfaffian version of the Binet-Cauchy theorem (see |3] for the details): 

Lemma 2.1. The T-functions of U.S]) with \l.f)]) can he expressed by 

Tn= Pf(S(n,...,i2n))Det(S(n,...,Z2„)). 

l<n<---<i2n<4A'' 

where E(ii, . . . , Z2n) is 2n x 2n suhmatrix of the 2n x AN matrix En, 

E{ii, . . . ,i2n) := 



Ei-^ ■ ■ ■ Ei, 



p,(2n-l) _ K;(2n-1) 



and B(ii, . . . , Z2n) is 2n x 2n skewsymmetric submatrix of the AN x AN matrix B, 



B{il, . . . , i2n) = 



Mr, 



V / 




The lower triangular part is given by bi^^i^, — —bi^^i^, and is left in blank. Also Det(i?) is given by 
the Wronskian determinant, 

Det{E{ii, • ■ • , i2n)) = WT{Ei^ E,^^^ ) = J| (k„ - pi^ ) exp 

m<n 

The sign is due to the order il.4\ ), *.e. pi < P2 < ■ ■ ■ < Pan ■ 

From this formula of the t function, the soliton solutions are completely determined by the nonzero 
coefficients of the pfaffians Pf(i3(«i, . . . ,i2n)) which are the Pliicker coordinates of the orthogonal 
Grassmannan Gr(2n, AN). For example, if all the coefficients are positive, then we have (2n, 4iV— 2n)- 
soliton solution in u{x,y,t), that is, we have 2n out-going line solitons in y — > c» and AN — 2n 
in-coming line solitons in y —>■ — oo (see As we can see from that in particular, if n = we 
have an 2A''-soliton solution similar to the KP equation, i.e. all 2A solitons are A-solitons, and all 
the interactions are of T-types. Then the main purpose of this paper is to identify soliton solution 
of the DKP equation with a certain combination of nonzero coefficients. 

Example 2.1. As the simplest example, the pfaffian Pi{B{ii, . . . ,^4)) is given by 

Pf (B(l, . . . , 4)) = 61,263,4 - 61,362,4 + 61,462,3 . 
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Note here that if the indices in each term have partial overlap, then the product takes the minus 
sign, e.g. (1,3) and (2,4). This is a key for the classification of four A-solitons (see Section^)). 

From Lemma [2. II one should note that each term in t„ in H1.2|l is given by the product of 2n 
exponentials Ek, i.e. 

2n n 

j=i k=i 

where {ij\j = l,...,2n}~{ik,jk\ik<jk^ fc = 1, . . . , n }, and they are the indices for the elements 

As we will show, this structure of the r-functions will be important to 



of the i3-matrix, i.e. 6. 



identify a soliton solution in the asymptotics y —>■ ±oo. 

2.1. The B- matrix. Since the coefficients of the r- function are determined by the S-matrix, we 
give some remarks on the structure of the matrix. For a generic AN x AN skew-symmetric matrix 
B, i.e. B G so(47V), there is a decomposition called skew-Borel decomposition |1|, 

B = LJqL'^ , 

where L Cz G, the group of invertible elements in the set of lower-triangular matrices with nonzero 
2x2 blocks proportional to identity along the diagonal, i.e. 



fai 


■•• 





\ 







ai ■ ■■ 


















2N 




* 






Y[a,^0 > 


* 


* 






i=l 


* 


* * 


* a2N 











* 


a2N/ 


> 



The matrix Jq is the 4A^ x 4A^ skew-symmetric matrix whose 2x2 diagonal blocks are given by 
Q j and all the other entries are zero, i.e. 



/o 



(2.1) 



1 

\0 -10/ 

A non-generic element of so(4iV) may be obtained by a permutation matrix tt G San with tt^^ ~ tt^ 
(i.e. TT G 0(4A^)), that is, we have 

B = ttLJqL^t:-^ = {T:L)jQ{7:Lf , 

where L is not a generic element in Q. 

Example 2.2. Consider the case = 1 (4 x 4 _B-matrix). Then the generic element B can be 
expressed hy B — LJqL'^ with 



L = 

Then the matrix B is given by 

(2.2) B = LJqL'^ = 









'\ 


a 








d 


/ 





e 





// 



eg. 



/O a2 ad 
~ab 


V 



ae 
—ac 
-cd + be + p 
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(We leave the lower triangular part in blank for skew symmetric matrix.) If a skew-symmetric 
matrix has zero at (l,2)-entry (i.e. 61^2 = 0) and others are nonzero, such matrix, say Bi, cannot 
be expressed in this form. In this case, one can consider the form B = {ttL)Jq{ttL)'^ with 





A 








0^ 




fa 








'\ 










1 










a 













1 










b 


d 


/ 















V 




\0 


e 





// 



The matrix Bi is then given by 



/O 



(2.3) 

In particular, with L 
(2.4) 



{■S2,3L)Jo{s2^3L)'^ 





ae ^ 


ab 


be + p 





—ad 




/ 



Id, the 4x4 identity matrix, we have Bi 

/O 1 0^ 
1 

Ji :— 



Ji defined by 




V 0/ 

Other non-generic element of _B-matrix, say B2, having 61^2 



bi,3 — can be expressed by B2 



(2.5) 



with 



Bo = 



/O 


V 








a 
ac 
—ab 

/ 



S2A 



(\ 






In particular, for L — Id we have B2 






-1 




1 



0/ 



and L — 



la 













a 








b 





/ 





V 











J2 defined by 



(2.6) 



J2 := 



/O l\ 
1 

V 0/ 

In the next section we show that those skew-symmetric matrices Jo, Ji and J2 define all the one 
D-soliton solutions. 

In general, one can start with Jq to define the r-functions, that is, r„ = Pf (Qn) with Q„ = EnJaSn 
of H1.6|) . Then consider a generalization of the Q„-matrix by replacing Jq with J^r := ttJott"^ for 
some TT g SiN n 0(4iV). The role of tt can be understood as the change of order of the column 
vectors in f„, i.e. 



(El, E2 



where the column vector E/c := {Ek,E^ 



, E4Ar)7r 

(1) 



(±E^(1), ±Ejr(2), 



,±E, 



(47V) 



) , and 7r(/c) indicates the permutation k 



Tr{k). The signs in w should be chosen properly so that the new r-function = (f„7r) Jo(£'„7r) is 
sign-definite, i.e. non-singular (see Section^ for more details). 

3. One D-soliton and two A-solitons 



Here we discuss the cases with 4x4 skew-symmetric matrix B, and show that the cases include 
one D-soliton and two A-solitons. 
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3.1. D-solitons. In Introduction, we show an example of one D-soliton H1.8|) whose r-function is 
given by ti = Pf(£iJo^j^) with Jq in (|2.1|l and £i in Lemma l2.ll for iV = 1. This D-sohton is 
denoted as [1,2 : 3,4], since it exchanges the asymptotic values of it; = (InTi)^: as (1,2) ^ (3,4) 
(recall (i,j) '■— Pi +Pj)- A general form of the one D-soliton can be expressed as [ii, ji : «2, ^2] where 
ii < 12 and ik < jk- There are three cases, 

[1,2:3,4], [1,3:2,4], [1,4:2,3], 

which correspond to the D-solitons generated respectively by the skew-symmetric matrices, Jq, Ji 
and J2 in Example 12. 21 as the B-matrices. Recall here that the soliton label [i,j : k,l] indicates the 
nonzero elements bij and bk,i in the _B-matrix. Namely for the J^-matrix with nonzero entries with 
1 at (ii, ji) and («2,j2) in the upper-triangular part, the r-function is given by 

(3.1) Ti = Qi^2 = {Pji -pi^)E,^Ej^ + (pj2 -Pi2)E^^Ej^ . 
Then the function w = (Inri)^, is 

9 1^1 1 

(3.2) w{x,y,t) := ^Inn = 2^^*^ + 2*-^'" ^ ' Pji)^^^^2^^^2,j2 -^*iji) 

fe=i 

with 9ij = 9i -\- 9j + In \pi ~ Pj\. For each asymptotic of x — > ±00, one of the exponential terms in 
the r-function (|3.1|l becomes dominant, and we have 



w{x,y,t) 



Pii + Pji 1 as a; — *■ —00 . 
Pi2 + Pj2 1 as a; ^ cxo . 



Thus one D-soliton exchanges (ii, ji) (*2, J2) with the values of w, and is denoted by [ii, ji : ^2, J2], 
as an element of the Weyl group of D-type. 

One soliton solution u = 2wx is a plane wave having the form, 

u{x, y, t) = (t){Kx + kyy - cot), 

Let us denote the wavenumber k = {kx,ky) and the frequency lo for (|3.2II by k[ii,ji : 12,32] and 
w[ii, ji : 12, j2], i.e. 

k[ji, ji : i2,i2] = {Pt2 +Pn ~ P^l " P31, Pf^ + Ph Pn ^P%)^ 

w[n,ii : ^2,^2] = -{pt_ +P% -Pn -P]i) 

The slope (or velocity) of the soliton in the x-y plane is given by c := dx/dy — —ky/k^ = —(j>f^ + 
Ph ~ Pii ~ P^) / iPi2 +-Pj2 ~ Pii ~Pji)- The peak of the soliton then determines the line in the x-y 
plane given by the equation Oi^.j^ (x, y, t) — Oi^j^ {x, y, t) for each t. 

We also note that the velocity c of one D-soliton can be obtained as follows: First note that 
a D-soliton of : k,l] gives a line 0ij = 6k,i- Then setting x = cy for this equation, we have 



u 



i,j - ^i,j{c)y + O^j with 



(3.3) "ij(c) ;= [Pt +P])c + pf +p] ■ 

Then taking the limit \y\ — > 00 for the equation 9ij = O^^i, the velocity c of the D-soliton is 
determined by Uij{c) = Uk^i{c), i.e. the intersection point of Uij and Uk^i- In Figure ^ we show 
Uij{c) of (|3.3|l for the parameter p :— {pi, . . . ,^4) = (—2, —1, 0, 3). Each solid dot gives the velocity 
of one D-soliton, and each open circle gives the velocity of an A-soliton, e.g. the point at ui^2 = ui^3 
gives the velocity of [2 : 3]-soliton. This implies that if 61^2 and &i_3 are only nonzero elements in 
the B-matrix, then ri gives one A-soliton of [2 : 3], i.e. n = (^1,2(^2 — Pi)E2 + ^1,3(^3 — _pi)£'3)£'i. 
Figure^can be also used to determine the solitons appearing from the ri-function in the asymptotics 
\y\ CO. For example, if bij > for all 1 < i < j < 4, then ri contains all six terms EiEj with 
i < j- Then from Figure ^ we have two A-soliton solutions of T-type, [1 : 3] and [2 : 4] (see below 
for the details). 
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' U12 

Figure 1. The functions Uij{c) and the velocities of sohtons. The intersection 
point between Uij and Uk,i gives the velocity of D-soliton : k,l]. The solid 
dots give the velocities of D-solitons, and the open circles give the velocities of A- 
solitons, e.g. the velocity of [1 : 2]-soliton is c = ^{Pi +P2) which is given by either 

Ul,4 = U2,A or Wi,3 = U2.3- 




[12:34] [2:4] [1:3] [2:4] [1:3] 




Figure 2. Resonant bifurcations of the D-soliton [1,2 : 3,4]. Each region, 
corresponds to a dominant exponential marked by the element bi_j. For example, 
with 62,3 7^ 0, [1,2 : 3, 4]-D-soliton bifurcates into [1 : 3] and [2 : 4] A-solitons as 
y ^ —00 (top figure). Those solitons can be found from Figure ^ 

3.2. Two A-solitons. One D-soliton can be considered as a degenerate case of two A-solitons (recall 
that A-soliton exchanges single letter as [i : j], while D-soliton exchanges two letters as [i,j : k,l]). 
Namely, D-soliton of : k,l] can be considered as the "product" of two A-solitons of [i : k] and 
[j : I] or [i : I] and [j : k], i.e. : k,l] = [i : k] ■ [j : I] — [i : I] ■ [j : k] as the product of two 

permutations. Then two A-soliton solutions can be obtained by adding extra exponential terms, i.e. 
extra nonzero entries to the _B-matrix. If one adds just one term, then we have Y-shaped solution 
satisfying the resonant condition. For example, for the _B-matrix having two nonzero entries bi-^j-^ 
and 6i2.j2 the upper triangular part, if we add a nonzero entry at either (ii,j2) or (12,^2), we 
obtain a resonant Y-shaped soliton with the resonant condition, 

f k[ii : i2] + k[ji : js] = HhJi ■ «2, ^2], 
\ uj[ii : 12] + uj[ii ■■ 32] = L^[h,ii ■ 12,32] ■ 

Here the wavenumber k[i : j] and the frequency Lo[i : j] are for the one soliton solution of the KP 
equation, i.e. 

k[ii : i2] ^ (p^i ~Pi2,p1^ -pD' ^[^i ■ ^2] ^ ~ P^, 
In Figure m (see also Example 13.11 below), we show the resonant bifurcation of [1,2 : 3,4]-soliton 
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into two A-solitons of [1 : 3] and [2:4]. Those solitons can be found from Figure ^ For example, 
in the left figure of Figure [3 we have 1/1.2,^2,3 and U3.4 from the nonzero entries of the _B-matrix. 
Then the three intersection points of the graphs of those Uij in Figure ^ give one D-soliton (solid 
dot) and two A-solitons (open circles). Figure ^ also shows that the solid dot indicates the pair of 
dominant exponentials (-E1.2, £"3,4) for [1,2 : 3,4]-soliton as y — > 00, while the open circles show the 
dominant pairs (i?i.2, £■2,3) for [1 : 3]-soliton and (i?2.3, £"3,4) for [2 : 4]-soliton a,s y ^ —00. 

Since there are at least four disconnected regions in the x-y plane divided by two line solitons, one 
needs to add at least two extra nonzero (positive) entries. Figure [21 shows the resonant bifurcations 
of a D-soliton into two A-solitons. Two soliton solution in this figure is a new T-type (referred to 
as TD-type) which does not exist in the KP equation (i.e. 7^ 0). The T-type of the KP equation 
can be obtained by a B-matrix with Det(i?) = 0, e.g. 

/O 1 2 l\ 
Oil 
1 

V 0; 

Note T2 — with this matrix due to Pf{B) — 0. One should also note that the generic skew- 
symmetric i?-matrix having all nonzero entries gives a two A-soliton solution of T-type (having a 
resonant hole), but it is not a solution of the KP equation. In this case, since those solitons are 
A-solitons, the function v^v^ = t^/'^i appears only locally in the x-y plane. Both r-functions for 
T-type two-soliton solutions of KP and DKP equations have six independent exponentials, i.e. Eij 
for 1 < i < j < 4, which is the key to produce T-type resonant interaction of two A-solitons (see 
also laiHl). 

Example 3.1. Consider [1,2 : 3, 4]-D-soliton which corresponds to the S-matrix having nonzero 
entries at (1,2) and (3,4), i.e. B ^ Jo in l|2.1|l : 

(a) [1 : 3]- and [2 : 4]-A-solitons (TD-type) by putting +1 at the entries (1,4) and (2,3). 

(b) [1 : 4]- and [2 : 3]-A-solitons (P-type) by putting -fl at the entries (1,3) and (2,4). 

Note that the TD-type in the case a) is not the same as the one of the KP-solitons, i.e. there is 
no resonance in this case. Original T-type of two A-solitons is obtained by the i?-matrix having 
all nonzero entries and Det(i3) = 0. Again there are three types of two A-solitons for DKP. These 
types play an important role for the classification problem as a building block. Figure |2| shows the 
diagrams of the three fundamental types of two A-solitons. The middle diagrams of 8-gon connect 
the types with the corresponding iJ-matrices as follows: We have nonzero bij if the pair (i, j) is 
connected in the diagram, and a corresponding A-soliton of [k : I] is given by the pair {k,l) having 
no connection in the diagram. Also note that the corresponding i3-matrix for 0-type is not in the 
generic one given in (|2.2|) , but in H2.3|) (compare this with the 0-type of KP solitons which belong 
to a Schubert cell of codimension one |Sj). Figure 01 shows those types of two A-solitons, which are 
the same as 2-soliton solutions of the KP equation (note here that we have set T2 — for all the 
types, i.e. ti gives a solution of the KP equation). The parameters are chosen as p = (—2, —1, 0, 3), 
and the velocities of solitons can be found from Figure ^ 



4. Solitons generated by 8 x 8 B-matrices 

Since each D-soliton has four phases, we need eight phases to describe a two D-soliton solution. 
The r2-function in the case with 8x8 _B-matrix is given by 



T2= J2 Pf(S(ii,...,z4))A(n,...,i4)exp I ^( 

l<n<-<i4<8 \k=l / 

where A(ii, . . . , 14) — Y[i<j<k<4iPij ~ Pik) > 0' ^^d the pfafhan is calculated as 

Pf (i?(il , . . . , ^4)) — ,*2 ~^ ,*4^*2 ,^3 " 
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1 O— O 2 

30— O 4 

0-type 





1 o- 



2 O — O 3 



P-type 



-O 4 



Figure 3. Three fundamental types of two A-solitons obtained from D-soliton. 
Each diagram provides the relation between i3-matrix and the corresponding two- 
soliton solution, i.e. bij ^ if and only if the pair is directly connected in 

the middle diagram, and [k : Z]-soliton implies that the pair (fc, I) has no direct 
connection. For example, in 0-type, two A-solitons are labeled [1 : 2] and [3 : 4], 
then we have the middle diagram to show the connections between the numbers 
in the pairs, i.e. (1,2) and (3,4) are not connected. This diagram implies we have 
nonzero elements &i,3, &i,4, 62,3 and 62,4. 






O-type 



~J\ 



Figure 4. Three fundamental types of two A-solitons. Those solutions have the 
same topological feature as the two soliton solutions of the KP equation (see also 
Figure 3 in [H|). The difference between those solutions appears as the existence of 
localized solutions (if u+w" — 0, then the DKP solution is also the KP solution). 
However note that the TD-typc in Figure 13 does not have a resonant quadrangle. 
The figures show the counter plots of the function w — (Inr)^; in the x-y plane 
for a fixed t, and the contrast of the shades shows the different values of w which 
increases in the positive x-direction. 



Let us now construct the _B-matrices for two D-soliton solutions which consist of only four nonzero 
entries in the upper triangular parts. We follow the steps: 

(1) Take four pairs {ik,jk), k — 1, .., 4 with ik < jk from the set {1, . . . , 8}, such that 

1 ^ ii < 12 < ■ ■ ■ < U, jk ji, [k ^ I). 

There are (2x4— 1)!! = 105 ways to make those pairs (compare this with (n~, n+) in the 
case of 4 A-solitons IS!). Those four pairs give four nonzero elements in the B-matrix, i.e. 
{ik,jk) = bi^j^ for fc = 1, . . . ,4. 

(2) Determine the signs of each pairs {ik,jk), so that all the pfaffians in the expansion of T2 
associated with the i3-matrix have the same sign. Since the pairs {ik,jk) are all different, 
each pfaffian has just one term, denoted by Pf((i,j), (fc, ?)) = ±bijbk,i- Here the sign depends 
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on the pairs, for examples, 

Pf((l, 3), (5, 7)) = 61,365.7, Pf((2, 6), (4, 8)) = -62.664,8- 

In order to determine the sign of the pfafRan Pf((i, j), (fc, I)) = zszbi jbi^ i, we define: 

Definition 4.1. We say that the pairs {ij,ik) and (ii,im) with ij < ii have a partial overlap, if 
ij < ii < ik < im- Otherwise, we say that the pairs have a won-partial overlap (i.e. total or no 
overlap). Let us introduce the sign of overlap between {ik,jk) and (ii,ji), 

Tkl - 



— if ik < ii < ik < ji (partial overlap) , 
+ otherwise . 

Then we have Pf((ifc, j^), {ii,ji)) = crkibi^.jkK,]i- 

Now we have: 

Lemma 4.1. Suppose that the signs Uki satisfiy 

0'l20'l3O'l4 = 0'l20'23O'24 — 0'i3Cr23O'34 — 0'140'240'34 • 

Then one can make the r-function T2 to be sign- definite, that is, all Pf(i3(ii, . 
sign. 



. ,44)) take the same 



Proof. We denote = sgn(ifc,jfc), so that sgn(Pf((ifc, jfc)(ii, j;))) = (Tu^k^i- The condition that 
all the terms have the same sign requires 

C^kl^k^l = CFk'V^k'^l' , 

for any k, I and fc', V . This leads to <Ji2<y3A = o'i3CT24 — ci4(T23 which gives the assertion. □ 
Applying the Lemma, one can prove: 

Proposition 4.1. Let i = (ii, ... ,14) and'} — {ji, . . . , J4). Then there exist 33 cases of sign- definite 
T-functions whose B-matrices have nonzero entries only at {ik,jk) o,nd {jk,ik) for k = 1,...,4. 
Those cases are given by: 



(a) 


For i — 


(1,2,3,4), 


there are 8 cases with j = (ji , . 


• • ,J4); 








(5,6,7,8), (5,8,7,6), (6,5,8,7), (6,7,8,5), 








(7,8,5,6), (7,6,5,8), (8,5,6,7), (8,7,6,5) 


(b) 


For i = 


(1,2,3,5), 


there are 4 cases, 










(4,7,6,8), (4,7,8,6), (6,7,4,f 


5), (8,7,4,6) 


(c) 


For i = 


(1,2,3,6), 


there are 4 cases. 










(4,5,8,7), (5,4,7,8), (7,4,5,f 


i), (8,5,4,7) 


(d) 


For i = 


(1,2,4,5), 


there are 4 cases. 










(3,8,6,7), (3,6,8,7), (6,3,7,f 


?), (8,3,7,6) 


(e) 


For i — 


(1,2,4,6), 


there is one case, (8, 3, 5, 7) . 




(f) 


For i — 


(1,2,3,7), 


there are two cases, (4,5,6,8), 


(6,5,4,8). 


(g) 


For i ~ 


(1,3,4,5), 


there are two cases, (2,6,7,8), 


(2,8,7,6). 


(h) 


For i — 


(1,2,4,7), 


there is one case, (6, 3, 5, 8) . 




(i) 


For i = 


(1,3,4,6), 


there is one case, (2, 8, 5, 7) . 




(j) 


For i — 


(1,2,5,6), 


there are two cases, (3,4,7,8), 


(4,3,8,7). 


(k) 


For i = 


(1,2,5,7), 


there is one case, (4,3,6,8) . 




(1) 


For i — 


(1,3,5,6), 


there is one case, (2,4,8,7) . 




(m) 


For i — 


(1,3,4,7), 


there is one case, (2, 6, 5, 8) . 




(n) 


For i = 


(1,3,5,7), 


there is one case, (2,4,6,8) . 





The case (n) corresponds to the i?- matrix B = Jq defined in H2.1|l . and the B-matrices for other 
cases can be expressed a.s B = Jj^ := ttJqtt'^ for some tt € Ss with appropriate signs in the entries of 
TT. Those signs are determined so that the r-function is non-singular. 
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4.1. Two and three D-solitons. We now classify D-solitons given by the r- functions associated 
with the S-matrices obtained in Proposition l4.1l We obtain: 

Theorem 4.2. Suppose we have the order, (ii,ji) < (i2ii2) < (^SiJa) < (H^ji) with (j-kijk) 
Pik + Pjk ■ Then as a generic situation, the D-soliton solutions of the DKP equation turns out to he 
the following two cases; 

(i) two D-solitons of [n, ji : 13,^3], [«2,j2 : «4,j4] or [ii,ji : 14^4], [i2,j2 : hJs], 

(ii) three D-solitons of[ik,jk ■ Ja], jz : ia,ja] and [imjm ■ ia,ja] for a = 2 or 3 (a,k,l,m 
are all distinct). 

Proof. First we note that each r2-function associated with the _B-matrix given in Proposition 
14.11 has six exponential terms. We denote them as EkEi for 1 < fc < Z < 4, where Ek '■= Ei^^j^ = 
exp((pij. + pji^)x + {pf^ + p'j^)y + ^fc) for a fixed t. To identify two dominant exponents which 
determine an asymptotic D-soliton, we set x — cy and define Ui,j :— {pi + pj)c + pf +p'j (see Section 
ISJ. Then for large \y\, the velocity c of asymptotic soliton can be given by the intersection point 
of 

'^ikjk — ^iijn where as a generic situation only two exponential terms become dominant, say 
EkEm and EiEm (see also for the details, and Figure [S] below). The soliton is then identified as 
[ikiik '■ ihji]- Then there are two cases where the number of solitons, i.e. the number of balancing 
pairs of dominant exponentials, is either two or three. To show this, we consider, for example, the 
case with i = (1,2,3,4), j = (5,6,7,8). We set the parameter p so that p^^^ < p'^^ < p''^^ < p^p 

fl") . . (2) (2) ("2") (2) 

whith p). :— pi^ +Pjk ■ Keeping this order, one can choose p, so that we have P4 < p^ < P2 Pi 
with p^^^ :— pf^ ~^ P'jk- (This can be done with a large negative value for pi and a small positive 
Ps.) Then it is easy to see from the graph of Uk ■= Ui^^j^, that one can have either two D-soliton of 
■■ hjs] and [12, j2 ■ 14,^4] or three D-soliton of [ii,ji : is, J3], [«2, J2 : and [13, J3 : 14, J4]- 

The parameters p = (pi, . . . ,P8) can take the same values except p^ for both cases (Figure gives 
an example for this, also see below). 

f2) f2) i1) f2l f2l 

One can also make a different order forp^ , e.g. p\ < P2 < Pi < P4 (by taking large positive 

Ps)- Then one can have three D-soliton of [ii, ji : i2, ^2], [*2, J2 '■ *3, js] and [12, J2 '■ *4,i4] with the 

velocities 02,4 < ci_2 < C2,3 where Ck,i is the velocity of [ikijk '■ *i, Ji]-soliton (recall that those are 

given by the intersection Uk = ui with Uk = Uik,jk)- 

(2) 

Thus there is a freedom in the ordering of , and this is a key to generate three D-solitons (com- 
pare with the case for the KP equation |21 IE) where we have only two A-solitons for six exponential 
terms in the r- function). 

In Figure |S1 we illustrate the case of i = (1,2,3,5) and j = (4,7,6,8): Here the ordering is 
(1,4) < (2,7) < (3,6) < (5,8). Then depending on the parameters p^'s, the number of asymptotic 
D-solitons is either two or three. The left (right) figure shows the case of two (three) D-solitons, i.e. 
there are two (three) values of c where two exponential terms become dominant. It is obvious that 
there is no case other than those two in the general case. □ 

We note here that all cases in Proposition 14. II have two- D-soliton solutions, but only some cases 
can have three-D-soliton solution. In particular, we have: 

Proposition 4.2. The case (n) in Proposition \4. 1\ can have only two-D-soliton solutions of [1,2 : 
5,6] and [3,4 : 7,8]. 

In order to prove the Proposition, we prepare the following Lemma: Let us first denote Ui^.j^ as 
Uk as before: 

ufe(c) := {p^, +Pj^)c + pf^ +p]^ , /c= 1,...,4, 
where {ik,jk) = (2k — 1, 2k) for the S-matrix of the case (n), i.e. B = Jq. Then we have: 

Lemma 4.2. Let Cj^k be the point of intersection given by Uj{c) ~ Uk{c). Then for a fixed index 
a E {1, ... ,4}, we have 

Ca.k < Ca,j < Ca,i , for i < j < k , a ^ i , j , k . 
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Figure 5. Two and three D-solitons associated with i = (1,2,3,5) and j = 
(4,7,6,8). The left figure shows two D-solitons with p = (-3,-2,1,2,3,4,5,6). 
The right one shows three D-solitons with the same parameters except ps — 8- The 
bottom figures show the graphs of Uij{c) = (jpi + pj)c + pf + p'j- The dots in the 
graphs show the dominant two exponential terms, which determine the asymptotic 
D-solitons. The squares correspond to the intermediate D-solitons forming the reso- 
nant quadrangle, and the triangles correspond to the resonant intermediate solitons 
from the pairs of three asymptotic solitons as shown in the top figures. 



Proof. From Mfe(c) = Uj{c), we have 

(P2j-1 + P2]) - {P2k-1 +P2k) 

Let us write P2k~i and p2k as 

P2k^l = Pk - Afc, P2k = Pk + Afc , 

that is, Pk = h{p2k-i +P2k), and we have 



c,,k = -{Pj + Pk) - 



Pj - Pk 

Suppose a < j < k (other cases follow the similar argument). Then we have 

A 2 _ A 2 /^2 _ ^2 

c.,~c.^k^Pk-P, + ^^-^^. 

It is then easy to show c„j — ^ > 0. □ 

Now we give a proof of Proposition 14.21 

Proof. From Lemma f4. 21 we have the order ci.2 < ci.3 < ci.4. Also note that the slope of Uk{c) 
is given by P2k~i + P2k —'■ Pk^\ hence p^^'' < p^2^ < pg^"* < p'"^\ Then we obtain qualitative graphs 
of Uk as shown in Figure |H| and from the graphs, we can find the dominant pairs of the exponential 
terms in T2 for y — » ±00. This shows that the dominant pairs are given by {EiE^^E^^E^) and 
{EiE4,EiE2) for y ^ 00, and (i?i£'2, £^2-^3) and {E2E3, E3E4) for y —00. This implies that we 
have two D-solitons of [1,2 : 5,6] and [3,4 ; 7,8] (recall Ek = E2k-i,2k)- This completes the proof. 
□ 

We also have: 
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Figure 6. The qualitative graphs of the functions Uk{c) = (p2fc-i +P2k)c+P2k-i + 
for the case (n) in Proposition 14.11 which states, e.g. ci_4 < ci^a < ci^2- The 
open circles correspond to two D-solitons of [1, 2 : 5, 6] with the velocity ci_3 and 
[3,4 : 7,8] with C2,4. 

Proposition 4.3. The case with i — (1,2,3,4) and j ~ (8,7,6,5) (one of the cases in (a) of 
Proposition \4. 1^ can have both two and three D-solitons. 

Proof. Follow the argument in the proof of Theorem l4.2l Also see the following examples: 

(i) For two D-sohtons, we take p {pi,p2, . . . ,ps) = (— 5, — 4, — 2, — 1, 0, 3, 3.5, 7). Then the 
open circles in the left figure in Figure [7| show the dominant exponentials which correspond 
to the D-solitons of [ii, ji : 14, J4] — [1,8 : 4,5] and [12,^2 : 13, js] = [2, 7 : 3,6]. 

(ii) For three D-solitons, we take p — (— 5, — 3.5, — 2, — 1, 0, 3, 4, 7). The three open circles in 
the right figure in Figure [7| correspond to the two D-solitons of [ii,ji : «2,j2] = [1,8 : 2,7], 
[«2,J2 : iajs] = [2,7 : 3,6] and [12, J2 : iaJa] = [2,7 : 4,5]. 

In FigurelHl we show the interaction patterns for those cases shown in Figure[7|(they are topologically 
the same as in Figure Isj. □ 

In Section O we will extend Propositions 14 . 21 and 14 . 31 to the general cases as Propositions 15 . II and 

Recall that a D-soliton can be identified as an element of the Weyl group of D-type, W^. In 
the case of two D-solitons, we have, for example, tti^s : (ii, ji) (^Sjjs) and 7r2,4 : (j2,j2) ^ 
(j4,j4). The group generated by those elements is an abelian subgroup of , and the orbit 
of this subgroup of («i, ji, 12, ^2) represents the asymptotic regions divided by the solitons, i.e. 
(«i,Ji,«2,j2), (ii,Ji,«4, j4), (12, j2,i3, js) and (z3,j3, 14,^4), which correspond to the dominant ex- 
ponentials in the r-function, that is, E1E2, E1E4, E2E3 and E3E4, respectively. In the r-function, 
we have two more exponential terms labeled by («i, Ji, is, js) and («2, ^2, «4, J4), i-c E1E3 a.nd E2E4. 
Those terms become dominant at the points of the resonant quadrangle as in the case of two A- 
solitons of T-type. Namely we have a resonant two D-soliton. The resonant condition is given by, 
for example, 

{Hh,ji ■■ iijs] +k[i3,j3 : 24,^4] = k[ii, ji : 24,^4] 
^^[hJi ■■ hJa] +^^[hj3 ■ iiji] = : 14,^4] 

In the case of three D-soliton, for example [ii, ji : 42, ^2], [«2,i2 : *37 ja] and [?2,i2 : *4,i4], we have 
7ri,2, 7r2,3 and 772,4 as the generating elements of the subgroup. Then the subgroup contains 7ri,3, 7ri,4 
and 773,4, and those are the D-solitons generated by the resonances of pairs, i.e. -Kkj represents 
D-soliton of [ik,jk ■ ihji]- Each resonant relation can be represented by the product 77^,2 ■T^2,j — T^ij- 

For any two D-soliton solution in Proposition l4.1l we have: 

Proposition 4.4. Two D-solitons constructed above are all in resonance, in the sense that there is 
a resonant quadrangle at the intersection point. Namely all the cases are of T-type with one hole. 
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Figure 7. The graphs of Ufc(c) = (p/c +P9_fe)c + p^, for fc = 1, . . . , 4. The in- 

dices are marked by {ik,jk) = (fc, 9 — fc) for fc = 1, . . . , 4. The left figure shows two D- 
sohton solutions of [1, 8 : 4, 5] and [2, 7 : 3, 6] with p = (—5, —4, —2, —1, 0, 3, 3.5, 7), 
and the right one shows three D-solitons of [1,8 : 2,7], [2,7 : 3,6] and [2,7 : 4,5] 
with p = (—5, —3.5, —2, —1, 0, 3, 4, 7). Note that only p2 is changed. 




Figure 8. Interaction patterns for two and three D-solitons in Figure [3 The 
function w = (lnr2)a; takes the values, w — > i2j J2) as x — * —00, and w —>■ 

(*3j jsj Hjjt) 81S X ^ 00. All six D-solitons in each figure can be identified as the six 
intersection points of the graphs Uk '■— Ui^j^ in Figure |7| 



Proof. First note that there are six exponential terms in the r-function. Each exponential term 
becomes dominant in some region in the x-y plane. For example, consider two D-solitons with 
[1,2 : 5,6] and [3,4 : 7,8]. Then the function w = (lnT2):r takes w — > (1,2,3,4) for x —00, 
and UI — > (5,6,7,8) for x 00. With those D-solitons, we have the regions marked by (1,2,7,8) 
and (3,4,5,6). In addition to those exponentials, we also have (1,2,5,6) and (3,4,7,8). As in the 
case of two solitons of the KP equation (see [2]), a resonant interaction of Y-shape can occur. For 
example, [1, 2 : 5, 6], [1, 2 : 7, 8] and [5, 6 : 7, 8] solitons are in resonant, and the three regions divided 
by those solitons are marked by (1, 2, 3, 4), (3, 4, 5, 6) and (3,4,7,8). Then the six exponentials in 
the T- function form two D-solitons in resonance (see the left figure in Figure IS]). □ 

Remark 4.3. In Figure 8 of [S], Isojima et al presented an example of (two) D-solitons which seems 
to have only one intermediate soliton (no resonant quadrangle). However, it turns out that this is 
not two D-solitons: First we note that this is the case (k) in Proposition 14.11 i.e. i — (1, 2, 5, 7) and 
j = (4, 3, 6, 8). This can be found by ordering the parameters (pi, . . . ,p4, gi, . . . , (74) in Figure 8 of 
[nj, i.e. relabel (pi, . . . , (74) in terms of our ordering (pi, . . . ,ps) withpi — > ps,P2 P7,P3 Pi, Pi — >• 
Pi,qi ^ pe,q2 ^ P5, 93 — > P3 and —^P2- The parameters chosen in the paper is then given by 
p = (— |, — |, — |, — i, i, |, |, 3). With those parameters, we note that pi + p^ =: (1, 4) — (2, 3) =: 
P2 +P3, a degeneracy in the parameters. Two D-solitons are given by [1, 4 : 5, 6] and [2, 3 : 7, 8], and 
the asymptotic values of w are (1, 2, 3, 4) and (5, 6, 7, 8) for x — > ^00. Then the resonant quadrangle 
consists of the D-solitons of [1, 4 : 2, 3], [1, 4 : 7, 8], [2, 3 : 5, 6] and [5, 6 : 7, 8]. However, because of the 
degeneracy (1,4) = (2,3), the [1,4 : 2,3]-soliton cannot exist, i.e. the function w cannot change the 
value across this soliton. Also notice that the [2, 3 : 5, 6]- and [1, 4 : 5, 6]-solitons are almost parallel, 
and so are the [1, 4 : 7, 8]- and [2, 3 : 7, 8]-solitons: For example, the peak of the soliton [2, 3:5,6] is 
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given by 6*2,3 = ^5,61 i-e- for a fixed t and with the values of p, we have the line for [2, 3 : 5, 6]-soliton, 

57 / 55 



(P5 +P6~P2~ P3.)x + + Pe - -P2 - P3)V V ^ 114^ j ^ 

which is almost parallel to the line 6*1 4 — 9c, q for the [1, 4 : 5, 6]-soliton, 

57 / 65 \ 

(P5 +Pe-Pi- Pi)x + [pI +pI-pI- pDv = ( ^ ^ 114^ j ^ 

Thus the soliton solution in Figure 8 of 5 consists of two D-solitons of [1,4 : 5,6] and [2,3 : 7,8] 
for y — » 00 and two other D-solitons of [2, 3 : 5, 6] and [1, 4 : 7, 8] for y —00, and the intermediate 
soliton is the [5,6 : 7, 8]-soliton. Other intermediate soliton of [1,4 : 2,3] may be considered to be 
located at y = —00. 

For any three D-soliton solutions in Proposition 14.11 we have: 

Proposition 4.5. Any pair of three D-solitons constructed above is in resonance to form an Y-shape 
vertex. This three D-soliton solution has no resonant hole and four resonant Y-shape vertices. 

Proof. Since there are six exponential terms which mark the six separated regions in the x-y plane 
bounded by those three solitons, this 3-soliton solution cannot have a resonant hole. Three solitons 
are labeled by, for example, [ii, ji : i2,.j2\, [*2, j2 ■ is, js] and [12,^2 : Hjji], and any pair of those 
solitons is in resonance, i.e. [ii,ji : «2, ^2] and [12, ^2 : *3, js] have a resonant interaction to generate 
[«i, ji; «3, is] -soliton (see the right figure in FiglHJ. It is easy to see that this resonance appear for 
any three solitons obtained in Proposition ^3 D 



4.2. Three soliton solutions consisting of two A-solitons and one D-soliton. Since one 
D-soliton is a degenerate two A-solitons, one can construct 3-soliton solution consisting of two A- 
solitons and one D-soliton from two D-solitons. One should note that resolution of one D-soliton 
into two A-solitons can work only for the case of two D-solitons. In the case of three D-solitons, the 
resolution of one D-soliton affects to other D-solitons, and it cannot produce two A-solitons. This 
can be checked by the asymptotic analysis using the graphs of Uij . 

Example 4.4. Consider two D-solitons corresponding to the _B-matrix of the case (n) in Proposition 
01 i.e. i = (1,3,5,7) and j = (2,4,6,8). Two D-solitons are [1,2 : 5,6] and [3,4 : 7,8]. Then, for 
example, two A-soliton of type [1 : 5] and [2 : 6] are obtained by adding two extra nonzero elements at 
(1, 6) and (2, 5). This fact can be verified easily from the diagram of T-type (up-side-down staircase 
with sides labeled by 1,2,5 and 6; see Figure|3J). The 3-solitons obtained here has only one hole (see 
Figure EJ). 

One can also get two A-solitons of [1 : 6] and [2 : 5] from [1,2:5, 6]-soliton by adding two nonzero 
elements at (1,5) and (2,6). Notice that those two A-solitons form the diagram of P-type (cannon 
shape with sides by 1,2,5 and 6; see Figure|2l. 

In the same way, one can resolve [3, 4 : 7, 8]-solitons into either the pair of [3 : 7] and [4 : 8] or the 
pair of [3 : 8] and [4 : 7] , by adding extra nonzero entries to the _B-matrix of the case (n) . 



4.3. Four A-soliton solutions from D-solitons. Those solitons can be constructed in the similar 
way as the case of three solitons (two A-solitons and one D-soliton) from two D-solitons. However, 
the compatibility among those A-solitons separated from D-solitons is somewhat complicated. For 
example, one can expect that two D-solitons, [1,2 : 5,6] and [3,4 : 7,8], generate two A-solitons 
[1 : 5],[2 : 6] from [1,2 : 5,6] and other two A-solitons [3 : 7] and [4 : 8] from [3,4 : 7,8]. The 
B-matrices for 3-solitons [1 : 5], [2 : 6], [3,4 : 7,8] and [1,2 : 5,6], [3 : 7], [4 : 8] may be respectively 
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12 56 
34 78 




Figure 9. Bifurcation of two D-solitons into three solitons with two A-soHtons and 
one D-sohton. The [1,2 : 5, 6]-D-soHton resolves into two A-solitons of [1 : 5] and 
[2:6]. D-soHtons are expressed as the double lines joined with the double circles. 
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Then adding 53,8,^4,7 to the _B-matrix in the left (or 61,6,^2,5 to the right one), the resulting B- 
matrix is expected to give four A-solitons [1 : 5], [2 : 6], [3 : 7] and [4 : 8]. However this matrix gives 
a singular solution, and a correct _B-matrix for this 4-soliton solution may be given by 



/O 



V 








1 

0/ 



Thus in this case we do not have a compatibility between the two sets of 3 solitons obtained from 
the D-solitons (in the sense that the r- function cannot stay to be sign-definite). A general rule to 
construct a i?-matrix for four A-solitons from two D-solitons seems to be complicated. Here we give 
a direct way to construct a i?-matrix for a given four A-solitons, based on the i?-matrices obtained 
in Proposition 14.11 

1) Give four A-solitons labeled by a pair of four numbers, i.e. 

^ {n^,...,n^), n+ = (n^", . . . , n|) 

where 1 = < ■ ■ ■ < and < . Then each soliton is labeled by [n^^ ; n'^] for 
fc = 1, . . . , 4 (the same as in the KP solitons 8 ). This implies that the asymptotic values of 
w — (InT)^, are given by X]fc=iPn* ^'^^ ^ ~^ ±00. 

2) Take a pair of two solitons, say [n~ : nf] and [nj : n'j'], and draw the corresponding 
diagrams shown in Figure 13 and do the same for the other pair. Here the choice of the pair 
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should be consistent with the asymptotic values. Then adjust the signs for each horizontal 
edges of the diagrams, so that all the pfaffians have the same sign. If all are satisfied, then 
the diagrams give the i3-matrix. 
Example: Consider four solitons with [1 : 2], [7 : 8] and any other two solitons with the indices 
from {3, 4, 5, 6}. Then one can construct the i?-matrix as follows: 

The diagram corresponding to [1 : 2], [7 : 8] is a backward C-shape (P-type shown in Figure |3Jl, 
having 1,2,7 and 8 on the vertical edge. This imphes we have nonzero elements for 61^7, 613, 62,7, 62,8- 
Now if the other pair is [3 : 6], [4 : 5], then the diagram associated to this is a cannon-shape (O- 
type) having 3,4,5,6 on the vertical edges. This gives the nonzero entries 63,4, ^3, 5, 64,61 ^5,6- In this 
example, all those entries in the upper triangular part can take +1 for a sign-definite r- function. 
Namely we have 

1 1\ 
1 1 
110 
10 
10 



0/ 

In the case of two solitons of T-type, [3 : 5], [4 : 6], one can use the B-matrix having the same entries 
except 63,5 = 2, and in addition 63 g = 1, 64 5 = 1. 

Notice that the nonzero entries has two groups corresponding to two D-solitons [1,7 : 2,8] and 
[3, 4 : 5, 6], and those groups have no partial overlaps. This structure is true for any cases satisfying 
the nonsigular condition, i.e. all the pfaffians have the same sign. (Recall that the i?-matrix with 
bi,8 = 62,7 = &3,4 = 65,6 = 1, i-e. i = (1,2,3, 5), j = (8,7,4,6), also gives three D-solitons of 
[1,'8:2,7]', [2,7:3,4] and [2,7:5,6].) 

For the same 4-solitons as the previous example, i.e. [1 : 2], [3 : 6], [4 : 5] and [7 : 8], one can also 
have the following i?-matrix which gives the same T2 function, 

1 1 \ 

110 

-1 -1 



0' 

1 1 



0/ 

This matrix can be constructed from the i?-matrix for i = (1,2,3,6) and j = (5,4,7,8), the case 
(c) in Proposition 14. II (i.e. this gives two D-solitons with [1,5 : 2,4] and [3,7 : 6,8]). In Figure [TUl 
we show the example of four A-solitons based on two B-matrices: The case a) corresponds to the 
S-matrix for i = (l,2,3,6),j = (5,4,7,8), and the case b) for i = (l,2,3,5),j = (8,7,4,6). We take 
p = (—4, —3, —2, —1, 1, 3, 4, 6) for both cases. 

As an example of four A-solitons, we consider the case with [1 : 7], [2 : 8], [3 : 5] and [4:6], i.e. 
there are two partial overlaps (two solitons of T-type). We first note that this set of 4-solitons can 
be split into two groups of 2-solitons having no partial overlaps, {[1 : 7], [2 : 8]} and {[3 : 5], [4 : 6]}. 
Since they are both TD-type, we have the TD-type diagram (up-side-down staircase). This implies 
we have nonzero entries 61.81 &i, 2, ^2, 7, 67,8 from the set {[1 : 7], [2 : 8]}, and 63,4, 63,67 64,5, 65,6 from 
the second set. However this i3-matrix does not give the 4-solitons. In Figure ITTl the left figure 
shows the solution given by this B-matrix, two D-solitons, [1,2 : 5,6] and [3,4 : 7,8] with some 
additional resonances at the intersection region of those solitons. This can be explained as follows: 
This i?-matrix can be considered as a deformation (bifurcation) from the B-matrix of two D-solitons, 
[1,2 : 5,6] and [3,4 : 7,8], and the r-function with this i3-matrix gives resonant D-solitons, [1,2 : 



(4.1) 



/O 




B 



(4.2) 



/O 




B = 
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a) b) c) 

Figure 10. Four A-solitons obtained from the _B-matrix of two D-solitons or three 
D-sohtons. Two D-sohtons in a) is given by the i?-niartrix with i = (1, 2, 3, 6) and 
j — (5, 4, 7, 8). Three D-solitons in b) is given by the i?-niatrix with i = (1, 2, 3, 5) 
and j = (8, 7, 4, 6). Both cases have the same p = (—4, —3, —2, —1, 1, 3, 4, 6). The 
fom- A-sohton solution in c) is given by the B-matrix 1)4. Note that the _B-matrix 
of 14. 2|) gives the same solution. 



3,4], [1,2 : 7,8], [3,4 : 5,6] and [5,6 : 7,8]. Then TD-type of 2-solitons [1 : 7], [2 : 8] is obtained 
from the bifurcation of [1, 2 : 7, 8]-soliton, which is one of the resonant solitons. Also note that the 
exponentials associated to (1,2,7,8) and (3,4,5,6) become dominant for two D-solitons for either 
y —>■ GO or y —>■ —oo. However those exponentials should not appear for this set of four A-solitons, 
that is, those four A-solitons cannot be obtained from this B-matrix. This can be fixed by taking 
following i?- matrix (i.e. change TD-type to T-type), 



/O 



(4.3) 



B 







V 



1 2 l\ 

1 1 

12 10 

110 

10 



1 



In Figure im the right fi gure shows the 4-solitons generated by this i?-matrix. 

We also mention that the r2-function with the generic 8x8 B-matrix having all nonzero entries 
gives four A-solitons of T-type, that is, they are the same type of the 4-soliton solution of the KP 
equation. Just notice that the number of exponential terms in T2 is (4) = 70, which is the same as 
the T- function for the KP equation (see [2])- Then the asymptotic analysis of the t- functions shows 
that four A-solitons are given by [1 : 5], [2 : 6], [3 : 7] and [4:8], which are of T-type. An example of 
the generic i?- matrix for nonsingular r-function is given by the skew-symmetric matrix having I's 
in all the entries in the upper triangular part. 



5. MULTI-SOLITON SOLUTIONS 

For the general case with arbitrary M = AN, we consider the _B-matrix with 2A'' numbers of 
nozero entries bi^j^. for fc = 1, . . . , TV and 1 = ii < • ■ • < zat, ik < jk in the upper triangular part. 
Although there is a scheme to find the i?-matrix for sign-definite rjv (i.e. having all the same sign 
in the pfaffian coefficients), the generalization of Proposition im seems to be difficult. However, we 
can show the following Propositions regarding to the number of D-solitons for two special cases: 

Proposition 5.1. The TM-function 1^1. ^) with AN x 4A^ B-matrix given by Jq of 1^2. 1\) can have 
only N D-solitons in the generic situation. 
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Figure 11. The 4-solutions in the right figure is obtained from the i?- matrices 
given by (|4.3|l . The left figure shows two D-sohtons obtained by the same matrix 
except bi^Y = 0; ^2,8 = 0, 63,5 = and 64^6 = (note the extra resonant structure at 
the intesection region of those D-sohtons. 



Proof. Apply Lemma l4.2l and find the dominant exponentials by following the argument of "level 
of intersection" introduced in [Sj (the main idea here is to find the dominant exponents from the 
graphs of Uij(c)). D 

The other case corresponds to the B-matrix, B = J^, where is defined by :— ttJqtt'^ with 
TT G San, the permutation given by tt = [2 : 4A^] • [4 : 4A^ - 2] • • • [2iV : 2iV + 2], i.e. 



/o 



(5.1) 



1\ 





1 








0/ 



Then we have: 



Proposition 5.2. The number of D-solitons generated by the Tj^-function J^l.'^) with AN x AN B- 
matrix given by of 1^5.1]) can be any number from N to 2N — 1, depending on the choice of the 
parameter p. 

Here instead of giving a precise proof, we demonstrate a case with M = 16 {N — A): In Figure 



1121 we show the graphs of Uk (c) — (pk + Pn- 



k)c + pl + 



for fc = 1, . . . , 8. With a proper choice 



of the parameter p = (pi, . . . jPie), one can find four different cases of intersection patterns of those 
lines. Then D-solitons can be found at the intersection points marked by open circles in Figure [T^ 
that is, the level of intersection of those points should be four (i.e. there are four lines above those 
points). For example, D-soliton [4, 13 : 5, 12] is obtained by the intersection of M4 and u^, and two 
dominant exponents for this soliton are EiE2E^Ei and EiE2E^E^ which are obtained from four 
lines above the intersection point U4 ~ (this is the definition of the level of intersection; see [2] 
for more details). We have the following four cases: 

(a) Four D-solitons (the top left figure) of [1, 16 : 8, 9], [2, 15 : 7, 10], [3, 14 : 6, 11] and [4, 13 : 5, 12] 
with the parameter, 

p = (-16, -14, -12.5, -11, -8.5, -6, -4, -2, 2, 5, 6, 9.5, 11.2, 13.5, 14, 17) . 

(b) Five D-solitons (the top right figure) of [1, 16 : 8, 9], [2, 15 : 4, 13], [3, 14 : 6, 11], [4, 13 : 5, 12] 
and [4, 13 : 7, 10] wfth 



p = (-16, -14, -12.5, -11.6, -8.5, -6, -3, -2, 2, 4, 6, 9.5, 12, 13.5, 14, 17) 
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Figure 12. The graphs of Ufc(c) = {pk + Pn-k)c + pj. +p\-j^k for fc = 1, . . . , 16. 
Each open circle of the intersections of two functions uu corresponds to a D-sohton 
given by T4. For eaxmple, the point for ui — us corresponds to the D-sohton of 
[1, 16 : 8,9]. The number of open circles depends on the choice of the parameter 
p = (pi, . . . ,pig). Notice that only and later uq are changed to increase the 
number of solitons (see the text for the explicit values of p). 

Note here that three D-solitons have the common index (4,13). The [2,15 : 7, 10]-soliton 
in the case (a) causes the resonant interaction with two D-solitons [2, 15 : 4, 13] and [4, 13 : 
7,10]. 

(c) Six D-solitons (the bottom left figure) of [1, 16 : 4 : 13], [2, 15 : 4, 13], [3, 14 : 6, 11], [4, 13 : 
5, 12] , [4, 13 : 7, 10] and [4, 13:8,9] with 

p = (-16, -14, -12.5, -11, -8.5, -6, -4, -2, 2, 5, 6, 9.5, 11.5, 13.5, 14, 17) . 

(d) Seven D-solitons (the bottom right figure) of [k, 17 — fc : 4, 13] for all fc = 1, . . . , 8, except 
fc = 4, and 

p = (-16, -14, -12.5, -11, -8.5, -7.5, -4, -2, 2, 5, 6, 9.5, 11.5, 13.5, 14, 17) . 

One should note that those p values do not give generic cases, since some of the lines Ufc(c) are 
parallel as shown in Figure [T^ However, one can make those to be generic by changing slightly the 
parameters p^'s without breaking the intersection patterns. (We show those non-generic parameters, 
since those provide a clear evidence of the existence of those solitons.) 

We also remark that all the A^-soliton solutions obtaind by ttv with the i3-matrix given in Propo- 
sition lOl are fully resonant cases like the T-type for the KP equation, that is, the r^r contains (^^) 

exponential terms having the form JlILi ~ Il^i ^ikJk (see 12] ). 

Finally, we note that in the case of generic i3-matrix having all nonzero entries, the rjv function 
has the maximal number of exponential terms {2^) which has the same structure of the r-function 
for fully resonant 2iV-soliton solutions (i.e. T-type) of the KP equation. Those 2N solitons are all 
A-type, and the interaction patterns are topologically the same for both KP and DKP equations. 
In those cases, the functions appear only locally at the interaction regions among 2iV-solitons. 
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